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Abstract 


The work described in this report details the application of extended 
Kalman filtering to estimating the Space Shuttle Propulsion performance, i.e. 
specific impulse, from flight data in a post-flight processing computer 
program. The flight data used includes inertial platform acceleration, SRB 
head pressure, SSME chamber pressure and flow rates, and ground based radar 
tracking data. The key feature in this application is the model used for the 
SRB's, which is a nominal or reference quasi-static internal ballistics model 
normalized to the propellant bum depth. Dynamic states of mass overboard and 
propellant burn depth are included in the filter model to account for real- 
time deviations from the reference model used. Aerodynamic, plume, wind and 
main engine uncertainties are also included for an integrated system model. 
Asstuning uncertainty within the propulsion system model and attempts to 
estimate its deviations represent a new application of parameter estimation 
for rocket powered vehicles. Illustrations from the results of applying this 
estimation approach to several missions show good quality propulsion 
estimates . 
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PRECEDING PAGE BLANK NOT FILMED 


WRITE ( 6 , 25) (MSID(J), J = I, ( I PACK ( 1 4 ) + 20 ) ) 


CALL ACCREQ ( 2 , I PACK ( 2 ) , DBNAME , I PACK ( 14 ) , MS ID , NUHBR , ISTAT , ERROR ) 
IF ( ERROR. NE.O ) THEN 
WRITE (6, 501) ERROR 

501 FORMAT ( ' ERROR IN ACCREQ - ERROR= ' , 15 ) 

STOP 
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80 CONTINUE 


WEHE ( M , L ) - (ALWS(M,L) + FCV ( M, L ) *DAR (M,L) ) * PX ( M , L ) * SQRT ( EM( L ) / 

(TX(M,L)*1545.0) ) * 3 . 8839 427 *< CD(H , L ) + F CV ( M , L) * DCD ( M , L ) ) 
WRITE ( * , 901) AREAHX ( L ) , RHO 
VEL a 14 4.0 *WEHE ( M f L ) / ( AREAHX ( L ) *RHO ) 
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GRAPH (I, 19) = VAL (31) 
GRAPH ( I , 20 ) = VAL( 32 ) 


DO 225 1=1, IMAX 

WHITE ( 9 , 995) SPTIME(I), (GRAPH(I,J), J = 1, 20) 
995 FORMAT ( 3(7E15.8,/) ) 

225 CONTINUE 
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YMAX = AMAX1 ( GRAPH ( I , NM ) , YMAX 
PVAR(I) = GRAPH! I, NM) 

250 CONTINUE 


CALL PLOT 2 ( XMIN, XMAX , YMIN, YMAX ) 

CALL PLOT3 ( TSHFT, PVAR(l), IMAX ) 

CALL PLOT4 < 9, 'TIME < SEC ) ' , 6 , YLAB , 1 3 , ' CONTRL INPUT 
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ACCESS ( IPACK , IW , MS X D , SRATE , TIME , TSTOP , STALE , VAL , I STAT 


CVAX SPALT * SNGL ( TIME 1/.3048 
SPALT = VAL(ll) 
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ROUTINE TO PLOT OUTPUT OF THE LFILTER.EXE PROGRAM 

DIMENSION YLBFIL(3 ,45) ,T<500) ,PVAR(500,3) ,XKM(500,30) ,PKM(500,30) 
1 ,RG(30) ,COVZ(500,15) , RESID(500, 15) 

DIMENSION YLAB(3) 

OPEN (UNIT— 2 , FILE”' YLABLFIL. DAT' ,STATUS-'OLD' , ACCESS- ' SEQUENTIAL ' , 
1 RECL-9 , CARRIAGECONTROL- ' LIST ' ) 

READ (2,5) ( (YLBFIL(I, J) ,1-1,3) ,J-1,45) 

FORMAT ( 2A4, A1 ) 


CLOSE ( UNIT-2 ) 

WRITE ( * , 6 ) ((YLBFIL(I,J) ,1-1,3) ,J-1,45) 
FORMAT ( T2 , 2A4 , A1 ) 


C 

10 

20 

C 

C 


C 


C 

30 

C 


C 

100 

C 


C 

C 

C 


C 


OPEN ( UNIT-1, FILE-' LFILOUT.DAT', STATUS-'OLD' , 

1 ACCESS— 'SEQUENTIAL' , FORM- ' UNFORMATTED ' ) 

OPEN ( UNIT-3, FILE- 'PLTFIL. DAT ' , STATUS-'NEW' ) 

I - 1 
CONTINUE 
J - 0 
CONTINUE 

J - J + 1 

READ ( 1 , END - 100) T(I) , TL, IMEAS , NMEAS 
WRITE ( * , 901) T(I) , TL, IMEAS, NMEAS 
READ ( 1 ) (XKM(I,K), K = 1, 30) 

READ ( 1 ) (PKM (I , K) , K - 1, 30) 

READ(l) (RG(K), K « 1, 30) 

READ ( 1 ) COVZ(I,J) 

READ ( 1 ) RESIDf I , J) 

IF ( IMEAS - NMEAS ) 20, 30, 30 
CONTINUE 

IMAX - I 
I - I + 1 

IF ( I.GE.500 ) GO TO 100 
GO TO 10 

CONTINUE 

XMIN - T ( 1 ) 

XMAX - T(IMAX) 

DO 200 KK - 1, 30 PRECLUftG FAQE BLANK NOT RLME'D 

DO 110 II - 1, IMAX 

SIGV = SQRT ( PKM(II,KK) ) 

BUPPER - XKM ( II , KK) + SIGV 
BLOWER - XKM(II,KK) - SIGV 

PVAR(II.l) - BUPPER 
PVAR(II,2) - XKM(II,KK) 

PVAR(II , 3) - BLOWER 
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original page is 
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c 

110 

c 

c 

c 


c 

c 


c 


120 

c 


c 

200 

C 

C 

c 

c 


c 

210 

C 

C 

C 


c 


c 


220 

C 


CONTINUE 

WRITE ( 3 , 903) ((PVAR(II. JJ) , JJ - 1, 3), II - 1, IMAX) 

WRITE ( * , 904) KK 
WRITE (*,6) YLBFIL(1,KK) 

YLAB(l) - YLBFIL( 1 , KK) 

YLAB(2) - YLBFIL(2,KK) 

YLAB(3) - YLBFIL(3,KK) 

WRITE (* , 6 ) (YLAB(L) ,L»1,3) 

YMIN - PVAR (1,3) 

YMAX » PVAR (1,1) 

DO 120 II - 1, IMAX 

YMIN - AMIN1 ( PVAR (II, 3) , YMIN ) 

YMAX - AMAX1 ( PVAR (11,1) YMAX ) 

CONTINUE 

CALL PLOT2 ( XMIN, XMAX, YMIN, YMAX ) 

CALL PLOT3 ( ' + ', T(l), PVAR(1,1), IMAX ) 

CALL PLOT3 ( T(l), PVAR(1,2), IMAX ) 

CALL PLOT3 ( T(l), PVAR(1,3), IMAX ) 

CALL PLOT4 ( 9, 'TIME (SEC) 6, YLAB(l) , 13 , 'PLTFIL ESTIMATE 
CONTINUE 

DO 300 JJ - 1, 15 

DO 210 II - 1, IMAX 

SIGR - SQRT { COVZ(II, JJ) ) 

PVAR (II , 1) - SIGR 

PVAR (II ,2) - RESID(II, JJ) 

PVAR (II , 3) - -SIGR 

CONTINUE 

WRITE ( 3 , 903) ( (PVAR (II , KK) , KK =1, 3), II = 1, IMAX) 

YLAB(l) - YLBFIL( 1 , JJ+30) 

YLAB(2) = YLBFIL(2 , JJ+30) 

YLAB(3) - YLBFIL(3, JJ+30) 

YMIN - PVAR (1,3) 

YMAX - PVAR (1,1) 

DO 220 II - 1, IMAX 

YMIN - AMIN1 ( PVAR (II, 3) , YMIN ) 

YMIN - AMIN1 ( PVAR (II, 2) , YMIN ) 

YMAX - AMAX1 ( PVAR(II.l), YMAX ) 

YMAX « AMAX1 ( PVAR(II,2), YMAX ) 

CONTINUE 

WRITE ( * , 904) JJ 

WRITE ( * , 6) YLBFIL(1, JJ+30) 

CALL PLOT2 ( XMIN, XMAX, YMIN, YMAX ) 

CALL PLOT3 ( T(l), PVAR(1,1), IMAX ) 
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c 

300 

C 

C 


C 

c 

901 

902 

903 

904 
C 


CALL PL0T3 ( ' + T(l), PVAR(1,2), IMAX ) 

CALL PLOT3 ( T(l) , PVAR(1,3), IMAX ) 

CALL PLOT4 { 9, 'TIME (SEC) 6 , YLAB( 1) , 13 , 'PLTFIL RESIDUAL 
CONTINUE 


CLOSE ( UNIT-1 ) 

CLOSE ( UNIT-3 ) 

STOP 

FORMAT ( 2E15.8, 215 ) 
FORMAT ( 5E15.8 ) 
FORMAT ( 5X, 3E15.8 ) 
FORMAT ( 2X, 15 ) 


) 


END 
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ARRAY! 1) = TIME 
BRRAY ( 1 ) = 0.0 
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IP13P9 = IP13 + 9 

IP22P9 = 1P13P9 + 9 
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310 CONTINUE 



DQ3 = X KM ( 9 ) - X ( 9 ) 

DQS = XKM ( 7 ) * * 2 + XKM<8)**2 + XKM(9)**2 
IF ( DQS. LT. 1.0 ) THEN 

DQ4 = — ( XKM( 7 ) *DQ1+XKM( 8 ) *DQ2+XKM( 9 ) *DQ3 )/SQRT ( 1.0 - DQS 
XKM (10) = X(10) + DQ4 
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FPLUME ( 1 ) = -( DFAB - FAORB + FAORB *FDP ) * FAB * RQ 

FPLUME ( 2 ) = 0.0 

FPLUME ( 3 ) = -DFNB*FABHDP 

TPLUME(l) = 0.0 



TPLUME ( 2 ) - DMMB * FABHDP 
TPLUME { 3 ) » 0 .0 
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TVEC ( 7 ) = — XKM( 7 ) /X (10) 

TVEC ( 8 ) = -XKM< 8 )/X( 10 ) 

TVEC ( 9 ) = -XKM< 9 )/X( 10 ) 

CALL MOLT ( TVEC, PKM, VEC, 1, NTS , NTS 
DO 170 I = 1, (NS+NPAR) 



320 CONTINUE 
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CALL NASSRB 

WRITE!*, 997) POHHAT ( IMTR ) 
WRITE!*, 997) YMSRB(IMTR) 


WRITE (2,997) YMSRB(II) 

CONTINUE 

GO TO 300 
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310 CONTINUE 
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SUBROUTINE SMLT ( C, A, B, NRA, NCA 
DIMENSION A(NRA,NCA) , B ( NRA , NCA ) 



END 

SUBROUTINE SYMTRK ( A, NS, AS 
DIMENSION A(NS,NS), AS ( NS , NS ) 
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-A ( 1 ) 



SUBROUTINE IMBED ( BIG, IM, IN, SUB, M, N, IROW, ICOL 
DIMENSION SUB ( 1 ) , BIG(l) 
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BIGA = A ( KK ) 



10 IF ( ABS(BIGA) - ABS ( A ( I J ) ) 
15 BIGA » A(IJ) 

L(K) = I 
M ( K ) = J 
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A ( I J J = H0LD*A(KJ) + A(IJ) 
65 CONTINUE 



A ( KK ) = 1.0/BIGA 
80 CONTINUE 
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DO 100 J = 1, 3 
AINV ( I , J ) = DET*B(I,J) 
100 CONTINUE 


END 

SUBROUTINE ZEROM< A, NROW , NCOL 
DIMENSION A(NROW,NCOL) 
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